Quadrature entanglement and photon-number correlations accompanied by 

phase-locking 
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We investigate quantum properties of phase-locked light beams generated in a nondegenerate op- 
tical parametric oscillator (NOPO) with an intracavity waveplate. This investigation continuous our 
previous analysis presented in Phys.Rev.A 69, 05814 (2004), and involves problems of continuous- 
variable quadrature entanglement in the spectral domain, photon-number correlations as well as the 
signatures of phase-locking in the Wigner function. We study the role of phase-localizing processes 
on the quantum correlation effects. The peculiarities of phase-locked NOPO in the self-pulsing insta- 
bility operational regime are also cleared up. The results are obtained in both the P-representation 
as a quantum-mechanical calculation in the framework of stochastic equations of motion, and also 
by using numerical simulation based on the method of quantum state diffusion. 
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I. INTRODUCTION 

In the presently very active field of continuous variable 
(CV) quantum information processing a challenging goal 
consists in generation of entangled intensive light beams. 
Various quantum optical schemes generating entangled 
bright light have been proposed for this goal. For the first 
time the CV entangled states of light have been studied 
in and demonstrated experimentally in for non- 
degenerate optical parametric oscillator (NOPO) below 
the threshold. An experimental progress in this direction 
has been recently reported in (Ü-El- However, up to now 
the generation of bright light beams with high level of 
CV entanglement meets scrious problems. One of these 
is degradation of entanglement due to uncontrolled dis- 
sipation and decoherence, and cavity induced feedback. 
Beside this, it is known, that the relative phase between 
subharmonics in above threshold NOPO undergoes dif- 
fusion process. This process destroys the frequeney de- 
generacy of modes and hence limits the production of 
quantum-twin beams in NOPO above threshold. For re- 
ducing such phase diffusion phenomena various methods, 
based on phase locking mechanisms have been proposed 

The simplest scheme proposed and experimentally re- 
alized by Mason and Wong [5j is the NOPO with addi- 
tional intracavity quarter- wave plate to provide polariza- 
tion mixing between two orthogonally polarized modes of 
the subharmonics. Indeed, it was shown that induced lin- 
ear coupling between generated modes results in a lock- 
ing phenomenon and hence in frequeney degenerate op- 
eration above threshold. Following this experiment, the 
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semiclassical theory of NOPO with self-phase locking was 
developed in Q and its further more detailed considera- 
tion has been presented by Fabre and colleagues [ïï| . In 
a more general approach the semiclassical theory of the 
self-phase-loeked NOPO has been investigated by Gro/3 
and Boller 

Type-II optical parametric oscillators with self-phase 
licking have recently attracted a lot of attention as efi- 
cient source of nonclassical light [ïl lïl lïl lïl lïl fl7| . 
A full quantum mechanical treatment of this system in 
application to generation of CV entangled states of light 
beams under mode phase-locked condition has been pre- 
sented in 0, 0| , where the regimes below, near and 
above threshold was considered. Quantum optical effects 
have been demonstrated by Fabre group in the series 
of experiments 0, 0, 0. Generation of quadrature 
EPR entanglement has been demonstrated experimen- 
tally with self-phase-locked NOPO below threshold. 
Experimental investigation of intensity quantum correla- 
tion with phase-locked NOPO above threshold has been 
also performed [ÏÍ|. The intensity-difference squeezing 
in electronically phase locked NOPO above threshold as 
well as the Hong-Ou-Mandel interferometry using twin 
beams have also been experimentally demonstrated by 
Feng and Pfister ■ Stable generation in application to 
metrology has been also demonstrated [ïsj . 

In this paper we continue investigation of phase-locked 
NOPO follow the paper 0. The previous studies of en- 
tanglement in the presence of phase localizing processes 
have only concentrated on stable stationary solutions for 
dynamics of intracavity subharmonics. Quite recently, 
the experimental observation of dynamical signatures of 
self-phase locking in a triply resonant degenerate OPO 
was reported in |19j for both stable and unstable regimes, 
however, the study of dynamical aspeets of entanglement 
for this system have been postponed for future investi- 
gations. In addition, here we will consider entanglement 
production in the unstable regime of generation. We shall 
also give below the quantum description of the unstable 
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regime in phase-space on the framework of Wigner func- 
tion. 

Our goal is two-fold. In one part of the present paper, 
we expand the previous study of phase-locked NOPO 
in the stable, steady-state regime of generation [Ï2| . 
particularly considering both the important details in 
phase-space in the frame of Wigner function and photon- 
numbers difference squeezing in the presence of phase 
locking. We also investigate quadrature EPR entangle- 
ment in the spectral domain in addition to our previous 
results which have been performed in the time domain. 
The other part of the paper is devoted to the problem of 
entanglement in the self-pulsing instability regime. 

The system under consideration is the NOPO with 
quarter wave plate inside the cavity, interacting with the 
thermal bath. Due to explicit presence of dissipation in 
this problem, one has to write the master equation for 
the reduced density matrix of the system, which within 
the framework of the rotating wave approximation and 
in the interaction picture is 



dp 
dt 



in 



i=l 



[2üipaf - afciip - paj 



(1) 



where 



H = hAidfai + iHE (aj — 03) 

i=l 

+ ihk(ü3afa2 — a^a\a 2 ) + hx{a^a 2 + aia 2 ), (2) 

where aj are the boson operators for the cavity modes 
lüí. The mode 03 at frequency u> is driven by an external 
field with amplitude E, while a% and a 2 describe subhar- 
monics of two orthogonal polarizations at the degenerate 
frequency w/2. The constant k determines the efficiency 
of the down-conversion process, while x describes the 
energy exchange between only the subharmonic modes 
due to the intracavity waveplate. It should be noted 
that maximally reached value for \ m the experiments 
is x/l — 0.21. However in the following calculations we 
consider x/7 in a more wide range keeping in mind future 
experimental achievements. We take into account the de- 
tunings of subharmonics Aj and the cavity damping rates 
7i and consider the case of high cavity losses for pump 
mode (73 3> 71, 72), when this mode can be adiabatically 
eliminated. However, in our analysis we take into account 
the pump depletion effects. We will solve Eq.Q in the 
framework of P-representation and stochastic variables 
on one side, and also by using the well known numeri- 
cal quantum state diffusion method (QSD), according to 
which, open quantum systems are represented by the en- 
semble of quantum trajectòries [2jj. The density matrix 
is restored in this case from the ensemble averaging over 
photonic Fock states. 

The paper is arranged as follows. In Sec. II we re- 
cali some theoretical points of self-phase locked NOPO 



and we present analysis of both photon-number quantum 
correlation and phase locking in the base of the Wigner 
functions of generated modes. Section III is devoted to an 
analysis of quantum fluctuations of both modes in below- 
threshold operational regime as well as to calculation of 
the squeezed quadrature variance. In Sec. IV we investi- 
gate the self-pulsing instability regime of NOPO on the 
base of quantum trajectòries and in the phase space. We 
also discuss there the CV entanglement in the unstable 
regime of generation. We summarize our results in Sec.V. 



II. SIGNATURES OF PHASE-LOCKING IN 
THE QUANTUM CORRELATIONS AND IN THE 
PHASE SPACE 

At first, we shortly discuss the phase- locking phenom- 
ena in the semiclassical theory of NOPO. The stochastic 
equations of self-locked NOPO for the complex c-number 
variables 014 and /3i corresponding to the operators ai and 
a~l , in the regime of adiabatic elimination of pump mode 
have the following form |ï^: 



dai . 
~dt = " (71 



+ (e — \a\a.2)P<2, — ix&2 + -Ri, (3) 



80! . 

-m = " (71 



tAi)/3i 

+ (e- \(3 l (3 2 )a 2 + ixP2 + R\. (4) 

Here e = kE/j^ and A = k 2 /js. The equations for a 2 
and 02 are obtained by changing the indexes (1) ^ (2). 
Rx, R 2 are Gaussian noise terms obeying the following 
correlations 



(R!(t)R 2 (t?)) = —{E - k ai a 2 )S(t - i'), (5) 
73 



(R\(t)RÍ(t')) = —(E - kp x p 2 )5(t - t'). (6) 
73 

It has been shown in 0, ^| that these equations in 
the semiclassical approximation have stationary stable 
solutions only if the following relation holds 



4x 2 AiA 2 > (71 A 2 - 72A1) 



(7) 



otherwise the system has only unstable solutions. 

It should be noted that the phase locking condition (Q) 
obtained here in the adiabatic approximation of triply 
resonant NOPO has a universal form and is realized for 
the various schemes of self-phase-locked NOPO. It has 
been obtained at first 0, 13 for the case of doubly reso- 
nant NOPO. As shown in pj this condition coincides with 
analogous condition obtained for triply resonant NOPO 
without adiabatic approximation. 
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Let us recali the results concerned to the phases and 
photon numbers for the stationary, stable regime of self 
locked NOPO for the cases of Ai = A 2 = A, and 
7i = 72 = 7 |l2j |. In above threshold regime e > Eth — 
y/(x~ l^l) 2 + 7 2 i the steady-state solution for the pho- 
ton numbers is: 



n = nio = n 2 o 



^/ e 2_( x -|A|) 2 - 7 



(8) 



for A > 0. For the opposite sign of the detuning, A < 
the mean photon numbers are given by the same Eq.JSJ), 
while the phases read as 



(pxo = ^-Arcsin- (\+ |A|) + tt ( k + \ 



<P20 = -Arcsin- (%+ |A|) + tt I k - - 



(10) 



while the phases have been found to be 



<Pia = f2Q = --Arcsin - (x + |A|) + irk, (9) 
2 e 



(k — 0,1,2,..). Analyzing the system with help of 
a linear treatment of quantum fluctuations in the P- 
representation we have arrived [l^ to the following cor- 
relators 



J 



5n + 

Sri- 
5tp- 



(Sn + ,S(p + ) 
(Sri-, StpJ) 



4n [7(7 + An ) + (x- | A |) 2 ] 
4An (7 + An ) V - 2A «o(x~ I A \)sign(A), 



1 ( 4n oX (x- I A |), 
4 | A | x V 2 Xlsign(A), 



-(7 



-An (x- 



A \)sign(A) 
-A7 



2x"/sign(A) 
2 1 A | (x- | A 



(11) 
(12) 



r 



where Sn± = Sn2 ± 8n\, 5ip± = 5if2 ± Sfi, and Srii (t) = 
rii (t) — riio and 8(fi (t) = ifi (t) — (piQ. These correlators is 
written in terms of the stochastic variables, namely be- 
tween photon-number sum and phase sum in the modes, 
as well as between photon-number difference and phase 
difference in the modes. What is now interesting for us is 
the dependence of both photon number quantum correla- 
tion and phase fluctuations on the wave plate parameter 
X- 



A. Photon-number quantum correlations in the 
presence of phase-locking 

Let us at first consider the quantum correlations in 
twin light beams generated in self-phase locked NOPO. 
Intensity correlations of twin light beams are usually 
characterized quantitatively by the quantum fluctuations 
of the intensity difference between the generated beams 
which is normalizcd to the corresponding shot noise level. 
In this way the intense twin beams quantum correlations 
have been experimentally observed several years ago in 
NOPO operated above its threshold (UH^. Up to now 
the several interesting applications of intensity correla- 
tion have been proposed and realized |23l |2J, [25| . In this 
scction we analyze shortly the numbcr-diffcrence squeez- 
ing in the presence of phase localizing process. Most 
published studies of quantum correlations have usually 
been performed in the spectral domain by measure of 
the corresponding squeezing spectra. However, here we 
restrict ourselves by calculation only the integral quan- 



tity. To characterize the photon-number correlation we 
address to the variance of the fluctuations in the photon 
number difference: 



R = ((ajdi 



4a 2 ) 2 ) 



(13) 



This variance is expressed through the stochastic vari- 
ables using the relationship between normally-ordered 
moments of the operators and the stochastic moments 
with respect to the P-function. In the balance case 
(fil) = (712) = (n) and in the framework of the stochastic 
variables we get 



R = 2(n) + (5n 2 _). 
Then, using the formulas l|12f> we obtain 



R = no 



!A|+x 



(14) 



(15) 



Thus, the variance normalized to the level of fluctuations 
for the coherent state is 



R R ^14- X ) 

RN -2^-2 {1 + jAr 



(16) 



For x — * 0, Rn = 1/2, i.e. the normalized variance 
reaches 50 percent relative to the quantum noise level 
of a coherent state, in agreement with the result of the 
paper [2(| devoted to an ordinary NOPO. As we see, 



for %/|A| -c 1 there is a small aggravation of the quan- 
tum correlation between the photon numbers of the cor- 
responding twin fields in a self-phase locked NOPO in 
comparison with the case of an ordinary NOPO. 

The effect arising from the photon correlation in each 
modes as well as between twin modes are usually investi- 
gated with the aid of the second-order correlation func- 
tions 

g tí = {afaf), g 12 = {a\a\aia 2 ), (17) 

where (i = 1,2). We calculate these correlations in the 
linear treatment of quantum fluctuations and in terms of 
the stochastic variables as 

,9ii = .922 = .9 = nl + l({6n%) + (Sn 2 _)), (18) 
gi2=nl + \{{5nl)-{5n*_)). 

It is easy to verify that the following relation between the 
correlation functions hold 



512 -.9- -(5nl)=g+^(l-j^), (19) 
while the correlator g equals to 

2 7 n , (X- |A|) 2 , n x , /on . 
■ 9 = " 0+ 4A (1 + 7 ( 7 + Ano) ) 'T (1 -R ) · (20) 

When the wavc platc is not included, the formula H19|) 
is transformed to the relation between the correlation 
functions of an ordinary NOPO [26| . 

As we see, the insertion of a polarization mixer slightly 
destroys the photon number correlations produced by an 
ordinary NOPO. Having discussed the influence of phase 
locking on photon- number quantum correlations, we now 
turn our attention to investigate the Wigner function of 
self-phase locked NOPO. 



B. Phase locking on the Wigner function 

It is clearly seen from the Eqs. (|l·l[) and Eas. 1)12(1 that 
for sufficiently small x the terms < 5n 2 + > , < 8n 2 _ > and 
< Síp 2 ^ > does not depend on x, while the term < Sipt_ > 
is inversely proportional to it. Indccd, from Eas. (|12jl we 
have 

(^_) = ^[7 2 -|A|(x-|A|)]. (21) 

Therefore, decreasing of x m the correlator leads to in- 
creasing of phase difference fluctuations, and at last with 
X = we arrive to well known fact of phase diffusion in 
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FIG. 1: The Wigner function of the self-phase locked NOPO 
in the stable regime for the parameters: A/7 = 0.1, A1/7 = 
A2/7 = 10, xll = 0-1) a ) e H = 5, (below threshold), b) 
e/7 = 10, (at threshold), c) e/7 = 11, (above threshold), d) 
e/7 = 11, but xll — 0-5. Averaging is over 3000 quantum 
trajectòries. 



NOPO without wave plate. This phenomenon also ex- 
hibits himself in the framework of Wigner function and 
we turn to the consideration of this point. 

It has becn found in from the general point of view 
that the subharmonic modes has the phase space symme- 
try properties: the Wigner functions W\ and W% of the 
modes have a two-fold symmetry under the rotation of 
the phase-space by angle ir around its origin, 



Wi (r, 9) = W l (r, 9 + n) . (22) 

Here r, 9 are the polar coordinates of the complex phase 
space. We show below the results of numerical sim- 
ulations in the framwork of QSD method. Figure 1 
plots the Wigner function of one on the subharmonics 
of self-phase locked NOPO over transient time and for 
all operational regimes. A clear self-phase locking in 
the transition through the generation threshold is seen. 
The Wigner function is almost Gaussian below thresh- 
old and is squeezed near the threshold. Above-threshold 
the Wigner function has two separated peaks due to the 
phase-localizing processes leading to the phase locking of 
the subharmonic modes. The distance between peaks in- 
creases with increasing the photon number (decreasing 
A). Note, that the decreasing of wave plate parameter x 
leads to less localization of these peaks (compare Fig.l(d) 
(X/7 = 0.5) with Fig.l(c) (x/7 = 0.1)). 
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III. QUADRATURE ENTANGLEMENT IN THE 
SPECTRAL DOMAIN 



6a.it) andft(í) = <5ft(í), (i = 1,2) 



In this section we continue investigation of CV entan- 
glement in the self-phase locked NOPO. In the previ- 
ous studies we have calculated the variances of the 
relevant distance V- — V (X% — X%) and the total mo- 
mentum V+ = V {Y% + ía) of the quadrature amplitudcs 
of two modes X k — -j= [a k exp (—i9 k ) + a k exp (i8k)] , 

Yk = 7f \- a k cxp (~ iG k) - a-k exp (i9 k )] ,(k = l, 2), where 
V(X) = (V 2 ) — (X) 2 is a denotation for the variance 
and 9 k is the phase of local oscillator for the k-th mode. 
The two quadratures X k and Y k are non commuting ob- 
servables. We have considered the variance in the time- 
domain and have quantificd the CV cntanglcmcnt by the 
inseparability criterion for the quantum state of two op- 
tical modes 01 



v = \<y + + v-)<i. 



(23) 



It has been demonstrated recently that time-dependent 
quadrature variance could be observed by means of time- 
resolved homodyne measurements |28| . This approach 
particularly allows for applications in time-resolved quan- 
tum information protocols. Nevertheless, so far squeez- 
ing as well as CV entanglement have been mainly demon- 
strated in the spectral domain and not in the time do- 
main. Therefore, we calculate here the variance V in 
the spectral domain and in a fully quantum mechanical 
treatment of self-phase locked NOPO below threshold, 
E < E t h- 

We use the stochastic equations (J3), (QJ linearized 
around the zero-amplitude solution = (3® = in 
the following matrix form for small fluctuations (Xi{t) = 



dL -> 
W = - FL + K 



(24) 



Here 



F 



A 



-B 



-B* A* 



(25) 



(5a 1 ,5a 2 ) T , 5(3 = (Sp l7 Sf3 2 ) T , R a 



where 5a 

(i?i,i?2) T , R/3 = {Rt- Rt) 1 are two-dimensional column 
vectors. The 4x4 matrix F is written in the block form 
with 2x2 matrices 



\ ix 7 + zA y ' 
The noise correlators are determined as 



(R a (t)Rl(t'))=D5{t-t') 
with the following diffusion matrix: 



1 

1 



D 



B 

È* 



(26) 



(27) 



(28) 



The correlation functions of the quantum fluctuations are 
obtained in the following form (we show results of ^1 
with changing some misprints) 



J 



5a (5a) 1 



5a (SP) 1 



2(S* 4 -4A 2 X 2 ) 

-2 



£ 



1 S 2 , -2 x aJ 



S 2 , -2 X A 



2(5 4 -4A 2 x 2 ) V- 2 ^ A ' S 



r 



x(s 2 - 

A(S 2 



2A 2 ), A(S 2 -2 X 2 
-2 X 2 ), x(S 2 -2A 2 



(29) 



where S 2 is introduced as 



S 2 



7 2 + X 2 + A 2 -e 2 - (30) 
These correlators are obtained in a steady state regime 
I 



and do not dependent on time. However, for spectrally 
resolved measurements of the cavity output fields we need 
in the correlation functions in the spectral domain. To 
calculate them we use the Fourier-transformed linearized 
equations for the stochastic amplitudes 



a{t) 



t a(uj)(L·, a(u>) 



/2tt 



'5(í)dí, 



(31) 



6 



/3(t) = 



/2tt 



'/3(í)dí, 



(32) 



which are obtained as given below 

«wóa(w) = -À6a(w) + ÈS(5(ui) + R a {uj), 
ilüS(3{lü) = B*óa{üj) - À*S(3(lü) + Rp{u). 



(33) 



Here, the nonzero correlation functions of noise terms in 
thc frequency space are 



(R(uj)R t (cu')) =D8{üj+J). 



We rewrite these equations in the more compact form for 

L{uj) = ( Jst) 1 as 



L(üj) = (F + iuI)- 1 R(uj) 



(35) 



and then using the correlators from Ea.í l34l0 . we explic- 
(34) itly calculate 



(L(lü)L t (lü')) = ((F + tujï)^ 1 R(lu)R t (lu')(F + ilú' Ï)~ 1T ) = 
(F + iwï)- 1 D(F + lüú' Ï)- 1T 8{uj + lü'). 

I 



(36) 



This result can be transformed to the following form As a result, the correlation matrices (|37[1 can be also 

written in the following block form 



(L(u)L J \u')) = (F + ujÏ)- 1 (F + uj , Ï)- 1 D5(uj + u)') = 
= {F 2 +uj 2 Ï)- 1 D5(üj + üj') (37) 

by using the formula 



(L{oj)L> (o/)> 



Go G\ 

G" * fi* 
i 



5{u + u'), (41) 



D(F + iüjI)- li = (F + ilúI)~ D. 



(38) where the matrices Go, G\ can be calculated from the 
equations 



The further simplification of this result is connected with 
the structures of the matrices F, D which are written in 
the block forms ffi jl .ffi fy . As a consequence we have 



F l + Iui 1 



al + bò —cò 
~c*ò a*ï + b*ò 



(39) 



where / is the identity matrix and the following notations 
are used 



a = (7 + íA) 2 - x 2 + £ 2 + u 2 , 
b = 2^x(7 + ïA), 
c = 27è. 



(40) 



J 





G^ 




f'ea \ 




GoJ 


1 = 1 


{0 eò ) 



al + bò —cò 



which can be transformed to the following form: 



,(42) 



aGo + bàGo — còG* = eò 
aGi + bòGt - còG* Q = 0. 



The solution of these equations is obtained as 



(43) 



Go 



:{{db* - ea*)ï + (da* - èb*)a),Gi = 
I 



d 2 - e 2 



(dl — eò), 



(44) 



where d — aa* + bb* — cc* and e = ab* + ba* 



The further calculation of the spectral variance using 
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VQr/2) 




2 3 4 5 



(à/y 



V&r/2) 





(b) 


\3 







2 3 4 5 



(à/y 




4 ■ ■ ■ ■ w/y 

2 4 6 8 10 

FIG. 2: The spectral variances versus w/7. The parameters 
are: a)A/7 = 0, x/7 = 0, 0.2, 0.5 (curve 1, 2, 3), e/e th = 0, 5. 

b) A/ 7 = 0, x/7 = 0,0.2,0.5(curve 1, 2, 3), e/e th = 0,8. 

c) A/ 7 = 0.2, 0.5, l(curve 1, 2, 3), x/7 = 0.05, e/e th = 0, 5 



the P-representation is Standard. A detailed description 
of the method can be found in [29|. We include also 
the output coupler transmissivity 27 for cavity output 
fields so that the frequency spectrum of the squeezing 
variance V(9, u>) corresponding to the integral variancc 
in the formula (|23|) is calculated to be 

V(6, u) = l + S (1) (6, u - A) + S {1) (6, -üj - A), (45) 
where 



S< 1 > = 2(Gi,ii(w) + G 1 , 2a (-w) 

7T 

+2Re[e™G*. 21 (üj)}), 



(46) 



and Gq 21 1 Gi,ix, G\^2 are the elements of the matrices 
Gg, Gi. The variance V(6,u) of the output quadrature 
amplitudes is symmetric in frequency around üj = 0, and 
equal to unity in the vacuum or coherent signal case, 
and can only reach its minimum value at the definite fre- 
quency. For the case of ordinary NOPO, if \ = 0, this 
result is coincided with the analogous one well known re- 
sult [2íJ . Typical results for self-locked NOPO are shown 
in Figs.2, for the various parameters x an d for 6 = tt/2 



chosen to minimize the noise level V(9, üj). The Standard 
quantum limit is shown in grey. We see that the minimal 
variance spectra remain less than unity for all freqüències 
only for the cases of small vàlues of the parameter x- 

Figures[2ïa),|2íb) plot the squeezing spectra in the ab- 
sence of the detunings (Ai = A2 = 0) for two vàlues of 
the pump field: e = 0.5s t h (Fig. Ufa)) and e — 0.8e t h 
(Fig. Hfb)), and for the various vàlues of the parame- 
ter x : X — (curves 1), x = 0.2 (curves 2), x = 0.5 
(curves 3). The excellent squeezing spectra centered at 
üj = 0, occurs near threshold for x = that is the case 
of an ordinary NOPO [23. We point out that there is a 
corresponding decrease in the squeezing as the parame- 
ter x increases. Note, that for phase-locked NOPO the 
minima of spectra are shifted form the zero frequency 
(see Fig. 2 (a), (b) curves 2 and 3) even in the case of 
zero detunings. Figures |2fc) plot squeezing spectra for 
the case of nonzero detuning. These results for x — 
are in a good agreement with the results obtained for an 
ordinary NOPO [Ü- 



IV. SELF-PULSING REGIME AND 
ENTANGLEMENT 

Now let us pay our attention to the case of classically 
nonstationary regime of generation, when the inequality 
@ does not vàlid. In this regime analytical treatment 
of semiclassical and quantum equations is complicated, 
therefore, we simulate Eqs. © and Q on one side, and 
use the QSD method to obtain numerical solution of Eq. 
||TJ on the other side. We find that the semiclassical solu- 
tion of the equation J3J , (@J , without the noise terms and 
for Pi — a*, exhibits the self-pulsing instability: the pho- 
ton number of intracavity modes oscillates periodically. 
We demonstrate this oscillations in Fig. 3 (a) (curve 1) 
for one of the modes. It appears that self-pulsing exists 
in the whole range of violation of the inequality (JJJ , and 
does not depend on the what parameter is changed to 
violate it. We present also the phase space trajectory 
of the semiclassical solution in Fig. 4 (a), which has the 
form of squeezed circle. The individual quantum tra- 
jectory of the single mode of self-locked NOPO for the 
same parameters as in Fig.3(a) is presented in Fig. 3 (b) 
as a quantum mechanical calculation on the base of QSD 
method. It is seen that trajectory repeats the oscillations 
of semiclassical solution. 

Due to the quantum noise, the averaged over ensemble 
of quantum trajectòries result does not exhibit oscilla- 
tions (see Fig.3(a), curve 2). The analogous situation 
takes place for various quantum systems, particularly, 
also for chaotic systems, where quantum trajectory re- 
flects the chaotic nonstationary behavior of the semiclas- 
sical trajectory, but ensemble averaged results has sta- 
tionary solution (see [ÜJ )• It has been shown in Ref . |3p| . 
that quantum chaos manifests itself, particularly, in the 
Wigner function: it repeats the shape of Poincaré sec- 
tion of semiclassical counterpart. So, it is expected that 
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FIG. 3: Time dependence of the photon number: (a) classical 
trajectory (curve 1) and quantum ensemble averaged result 
(curve 2), and (b) quantum trajectory of self-phase locked 
NOPO in the regime of self-pulsing. The parameters are: 
A/ 7 = 0.1, Ai/ 7 = 10, A2/7 = -5, x/7 = 0.5, e/7 = 4. 
Curve 2 on the Fig.3(a) involve averaging over 12000 quantum 
trajectòries. 
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FIG. 4: Semiclassical phase space trajectory (a) and contour- 
plots of the Wigner function (b) of the self-phase locked 
NOPO in the regime of self-pulsing. The parameters are: 
A/7 = 0.1, A1/7 = 0.1, A2/7 = -0.1, x/7 = 0.5, e/7 = 3. 
Averaging in Fig.4 (b)is over 3000 quantum trajectòries. 



the features of the phase trajectory (Fig.4(a)) will be re- 
flected in the Wigner function also for the system under 
consideration. Indeed, as it is seen from Fig.4(a) and 
Fig.4(b), where the contour-plots of the Wigner function 
is depicted, it is really the case. Note, that the result of 
Fig.4(b) is obtained in the framework of QSD method. 
Thus, we can conclude that the photon number of self 
locked NOPO is stationary in the full quantum treat- 
ment, even when the semiclassical counterpart in non- 
stationary. However, self pulsing manifests itself in the 
Wigner function, which is rather diffcrcnt from the case 
of stationary stable regime (see Fig.l and Fig. 4(b)). 

It is also interesting to compare the result of Fig.4(b), 
which has a slightly squeezed circle form, with the 
Wigner function of ordinary NOPO (\ = 0), presented in 
the paper [3]]], which has the form of circle. This form is 
the manifestation of phase diffusion. Indeed, the phase in 
phase space is proportional to the angle between the Ox 
axe and radius- vector of the point. As the Wigner func- 
tion distributed uniformly to all directions, this mean 
that the phase operator is cntirely uncertain, which is 
equivalent to phase diffusion. For the case of self locked 
NOPO the distribution is almost uniform to all direc- 
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FIG. 5: Entanglement versus pump amplitude. The parame- 
ters are: A/7 = 0.1, A1/7 = 10, A2/7 = -10, \h = 0.1. 



tions, which means that the modified phase diffusion also 
takes place here. The squeezed form of the circle is due 
to oscillation of photon number. 

Analogous results, but for triply resonant cavity have 
been obtained numerically by Gro/3 and Boller pfl in the 
framework of semiclassical treatment: the phases of the 
modes of self locked NOPO in the self-pulsing regime 
diverge. Moreover, the rate of the phase diffusion is the 
same as for an ordinary NOPO. It is interesting that 
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FIG. 6: Entanglement versus wave plate parameter \H 
(squares - mimerical result, line - quadratic fit). The pa- 
rameters are: A/7 = 0.1, A1/7 = 10, A2/7 = —10, e/7 = 3. 

the semiclassical phase space trajectory in Fig.4(a) has 
similar form as in the case of triply resonant NOPO. 

Another interesting question, which arises from study- 
ing nonstationary regime is: whether entanglement takes 
place in this regime or not? And, if yes, does it reach 
the same extent or perhaps is higher with respect to the 
stationary regime. Nevertheless, we understand that it 
will not have the desirable property of sclf-phase lock- 
ing. For answering these qüestions, we consider the de- 
pendence of dispersion V H23J) on both driving field am- 
plitudc in units of its threshold value e/ethi ancl quar- 
ter wave plate parameter x- A.t first, we numerically 
estimate that £th/l = 1 f° r the parameters used, and 
then we consider the dependence of V on e/eth- Then, 
we perform calculations and the results are presented in 
Fig.5. It is seen, that the variance drops down near the 
threshold and achieve the minimal value V — 0.55 at the 
threshold. Then, it increases by increasing the driving 
amplitudc. This behavior is similar to that for the case 
of stable operational regime of self-phase locked NOPO 

m 

We also present here the dependence of V versus wave 
plate parameter \ in Fig. 6, where the decreasing of en- 
tanglement degree by increasing of \ is clearly evident. 
This means, that while, as mentioned in Sec. II, the local- 



ization of phases of the modes is improved with increasing 
of x, the entanglement is worsening. 

It is remarkable, that this dependence is exactly 
quadratic: we fit the numerical results (squares in Fig. 6) 
with quadratic curve and find excellent coincidence (sòlid 
curve in Fig. 6). 



V. CONCLUSION 

In conclusion, we have studied in the full quantum me- 
chanical manner the properties of the light beams gen- 
erated in self-phase locked NOPO. We have continucd 
the recent investigations of the quantum aspeets of this 
device in the steady state, stable regime of generation 
[Ï3 | in one side and also we have considered the specific 
signatures of the NOPO containing a birefringent ele- 
ment in the unstable regime of generation on the other 
side. In this latter case, the system exhibits the both 
sclf-pulsing temporal behavior and a new type of phase 
diffusion as it has been demonstrated by numerical sim- 
ulations on the framework of the Wigner function. The 
entanglement does oceur in this regime too, but suffcr 
some worsening: it decreases quadratically when wave 
plate parameter increases linearly. Considering the sta- 
ble regime of generation, we conclude that significant 
noise reduction in intensity difference variance as well 
as in the quadrature squeezing spectra are possible for 
small parameters x/l « !■ A detailed analysis of two- 
mode squeezing spectra for self-phase locked NOPO be- 
low threshold has been given, using a P-presentation in 
applications to the experiment recently performed |l4j . 
We have also analyzed photon-number correlations in the 
presence of phase localization which has also been stud- 
ied on framework of the Wigner function. 
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